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Abstract 

The spinless Salpeter equation is the combination of relativistic kinematics 
with some static interaction potential. The nonlocal nature of the Hamiltonian 
resulting from this approximation renders difficult to obtain rigorous analytic 
statements on resulting solutions. In view of this unsatisfactory state of affairs, 
we calculate analytic upper bounds on the involved energy levels, and, for the 
Coulomb potential, the ground-state energy at the critical coupling constant. 
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1 Introduction: The Spinless Salpeter Equation 

The most straightforward generalization of standard nonrelativistic quantum theory 
towards reconciliation with requirements imposed by special relativity is to describe 
a quantum system under consideration by the well-known spinless Salpeter equation. 
Consider a quantum system the dynamics of which is governed by a — by assumption 
self-adjoint — Hamiltonian H of the form 

H = T + V, (1) 

where T denotes the square-root operator of the relativistic expression for the kinetic 
energy of some particle of mass m and momentum p, that is, T = y/p' 2 + m 2 , and 
V = V(x) represents an arbitrary coordinate-dependent static interaction potential. 
The eigenvalue equation for this Hamiltonian H, H\xk) — Ek\Xk)i k — 0, 1, 2, . . ., for 
Hilbert-space eigenvectors \xk) corresponding to energy eigenvalues 

E = (Xk\H\xk) 
(Xk\Xk) 

is the one-particle spinless Salpeter equation. However, because of the nonlocality of 
H, it is extremely hard to obtain analytic statements from this equation of motion. 

Of particular importance in physics are "central" potentials, i.e., potentials which 
depend only on the radial coordinate r = |x|. Among these, the most prominent one 
is the Coulomb potential Vc(r), which is parametrized by some coupling constant a: 

V(x) = Vb(r) = -- , a > . (2) 

The semirelativistic Hamiltonian ([!]) with the Coulomb interaction potential V in (0) 
defines the spinless relativistic Coulomb problem. In the past, the interest in this has 
undergone an eventful history. (For a rather comprehensive review, consult Ref. Jl|.) 

By examining || the spectral properties of the Hamiltonian ([I]) with the potential 
(0) one infers the existence of a Friedrichs extension up to the critical value a c = 2/7T 
of a. As far as analytic statements about this spinless relativistic Coulomb problem, 
in particular, about the resulting energy levels, are concerned, up to now one has to 
be content with few series expansions of the energy eigenvalues in powers of a |3|, |J , 
which then are, of course, only significant for some region of rather small values of a. 



2 Variational Upper Bounds on Energy Levels 

Variational upper bounds on the eigenvalues of the above Hamiltonian H are derived 
rather easily in the following way ||. We introduce an arbitrary real parameter [i and 
make use of the positivity of the square of the (obviously self-adjoint) operator T—/jl, 

< (T — /i) 2 = T 2 + /i 2 -2fiT = p 2 + m 2 + /i 2 - 2/xT , 

in order to find, for the kinetic energy T, the operator inequality (see also Ref. M) 

T< p2 + m2 + ^ forah>>0, 
2 fx 

and, therefore, for the semirelativistic Hamiltonian H in ([I]), the operator inequality 

H < H s (ji) for all fi > , 
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with some Schrodinger-like Hamiltonian, denoted by H s (fi) and given by the unusual 
form 

2 n 

Now, invoking the min-max principle, we infer that the set of energy eigenvalues Ek, 
fc = 0, 1, 2, . . ., of the Hamiltonian H, if ordered according to E < E 1 < E 2 < . . ., is 
bounded from above by the corresponding set of energy eigenvalues Es } k(fj) of Hs(f-i), 
if the latter are similarly ordered according to £7s,o(a0 < Es,i(f-i) < Es^i^) < • • •, i-e., 
Ek < Es^ifi) for all fi > 0, and is, consequently, also bounded by the minimum of all 
these energy eigenvalues: 

E k < mm. E S k (fjL) . 

At>0 

For the Coulomb potential (§), these energy eigenvalues Es jn ([i) are given by 

1 



2 /i 



m 2 + /i 2 1 - 



a 2 * 



n 2 



where the total quantum number n is defined according to n — n t +£+l. Minimizing 
the latter expression with respect to the parameter //, we obtain as minimal bound H] 



/ a 2 

mm£ s , n (//)=myi-— . (3) 

Comparing with the results obtained in the framework of perturbation theory [|J , 
which read, for instance, for the level n = 2, I = 1, 

E _ _ a 2 _ 7 a 4 727 a 6 

m~ ~ T ~ 384 + 82944 + ' ' ' ' 

we observe that our upper bound (|3|) not only possesses the elegance of a particularly 
simple form but is also valid for all a < a c and arbitrary levels of excitation. A quick 
glance at Table [l] shows that the relative error of our bound is always less than 0.1 %. 

Table 1: Comparison of the perturbatively computed [|J energy eigenvalues for, e.g., 
the level n = 2 and £ = 1 with our variational upper bounds, extracted from Eq. (|3p 



a 


Perturbation Theory [| 


E/m 

Variational Procedure, Eq. (|3]) 


0.0155522 


0.999969765 


0.999969766 


0.1425460 


0.997452 


0.997457 


0.2599358 


0.99147 


0.99152 


0.3566678 


0.9838 


0.9840 


0.4359255 


0.975 


0.976 


0.5 


0.967 


0.9682 
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3 The Energy at the Critical Coupling Constant 

Let us apply again the min-max principle. Our particular (somewhat sophisticated) 
choice for the basis vectors \4>k) to be adopted here is defined by trial functions ipk{f) 
given in configuration-space representation (with the radial coordinate r = |x|) by @ 



■r fc+ P _1 exp(-mr) , P>0, m > 



\ 47rr(2fc + 2/3+ 1) 



Whereas k indicates a positive integer, k = 0, 1, 2, . . ., the parameter /3 is introduced 
to allow, for some given value of the coupling constant a, of a total cancellation of the 
divergent contributions to the expectation values of kinetic energy T and interaction 
potential Vc(r): (3 = /3(a). More precisely: in order to provide for that cancellation, 
the parameter /3 has to be adjusted for the ground state according to the relation || 

/7T 

a = p cot I — p 

which implicitly determines (3 as a function of the coupling constant a. In particular, 
this relation tells us that the critical coupling constant, a c , is approached for j3 — > 0. 
From the behaviour of the matrix elements (ipi\T\i/jj) of the kinetic energy T for large 
momenta p and of the matrix elements (ipi\Vc(r)\ipj) of the Coulomb potential Vc(r) 
at small distances r, respectively, it should become evident that, for our particular 
choice of basis states \ipk), these singularities can arise only in matrix elements taken 
with respect to the ground state \ipo), that is, merely in (■0 o |T'|V'o) an d {i>o\Vc(r)\xjjo) . 

The remainder of our way is straightforward. With the above basis functions, we 
obtain, by calculating the 2x2 energy matrix (and performing very carefully the limit 
(3 — > 0), for the roots of the corresponding characteristic equation the expression || 

E 2 



m 15 Ti 



60 In 2 -23 ± J (60 In 2) 2 - 4800 In 2 + 1649 



which yields the upper bound E /m = 0.484288 . . . for the ground-state eigenvalue of 
the Coulombic Hamiltonian. In principle, the d roots of the characteristic equation 
may be determined algebraically up to and including the case d — 4. For d — 4, our 
method constrains the lowest-lying energy level by the bound Eo/m = 0.4842564 . . . 
On the other hand, at the critical coupling constant a c the best numerically obtained 
bounds are given by || 

E 

0.4825 < — < 0.4842910 for a = a c . 

m 

Comparing these only numerically determined bounds on the ground-state energy Eq 
with the above upper bound, we realize that even our 2x2 analytical bound lies well 
within the numerically obtained range. Consequently, we obtain a clear improvement 
of the best upper bound available up to now for the ground-state energy level of the 
spinless relativistic Coulomb problem at the critical value of the coupling constant. 
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